Electric field dependent structural and vibrational properties of the Si(100)-H(2xl) 
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We report a first principles study of the structure and the vibrational properties of the Si(100)- 
H(2xf) surface in an electric field. The calculated vibrational parameters are used to model the 
vibrational modes in the presence of the electric field corresponding to a realistic STM tip-surface 
geometry. We find that local one-phonon excitations have short lifetimes (f ps at room temperature) 
due to incoherent lateral diffusion, while diffusion of local multi-phonon excitations are suppressed 
due to anharmonic frequency shifts and have much longer lifetimes (fO ns at room temperature). 
We calculate the implications for current induced desorption of H using a recently developed first 
principles model of electron inelastic scattering. The calculations show that inelastic scattering 
events with energy transfer nhoj, where n > 1, play an important role in the desorption process. 
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I. INTRODUCTION 

STM induced desorption of hydrogen(H) from the 
monohydride Si(100) surface offers the possibility of 
lithography with atomic resolution [y-^]. Investigations 
of the desorption mechanism ^[T^] have established the 
dependence of the desorption rate on the bias voltage, 
tunnel current and H isotope. At high positive biases, 
Vb > 4 V, the experimental results are consistent with 
electron induced desorption due to direct excitation of 
the Si-H bond by a single electron At negative 

and low positive biases the desorption rates show power- 
law dependencies on the electron current consistent 
with a multi-electron process j|,|l4|,[l5| , and the measured 
desorption rates are in quantitative agreement with first- 
principle calculations |ll|,0]. 

The desorption by multi-electron scattering is only 
possible because the H stretch frequency has a long life- 
time. The long lifetime is a result of a vibrational quan- 
tum too low for coupling with electron-hole excitations, 
while well above the Si phonon spectrum, and can thus 
only transfer energy to the substrate via a multi-phonon 
process. At room temperature experimental estimates of 
the lifetime due to this process are r w 10 ns pqjLOl , 
However, a local excitation is not an eigenmode and will 
decay into a H surface phonon by a coherent process. 
This decay is several orders of magnitude faster than 
multi-phonon energy relaxation, and must therefore be 
included in the theoretical models. It has been proposed 
by Persson and Avouris that the vibrational Stark 

shifts due to the electric field from the tip can localize 
vibrational modes below the tip. The localized modes 
may still transfer energy laterally by incoherent diffusion 
(the Forster mechanism) but it was found that this decay 



channel is also reduced by the Stark shifts. However, the 
work assumed that the electric field from the STM tip is 
localized on a single H atom below the tip, and this is 
not the case for realistic tip geometries. 

In this work we present results for the vibrational prop- 
erties of the Si(100)-H(2x 1) surface in the presence of the 
electric field from a more realistic model for the STM tip. 
The tip is described by sphere of radius, Rt — 500 A, with 
a protrusion of atomic dimensions, and we determine the 
electric field by solving the Poisson's equation numeri- 
cally. To obtain the effect on the H vibrations we set up 
a phonon Hamiltonian with parameters obtained from 
a first principles calculation of the vibrational proper- 
ties of the Si(100)-H(2x 1) surface in an external electric 
field. We find that the electric field does give rise to 
a localized vibrational state below the tip; however, its 
lifetime is very short (10 ps) due to incoherent exciton 
motion. However, we find that the anharmonicity of the 
Si-H bond potential reduces the lateral energy transfer of 
higher excited excitations(n > 1) of the Si-H bond. We 
present first principle calculations of the desorption rate 
taking this effect into account, and find that two phonon 
excitations play an important role in the desorption pro- 
cess. 

The organization of the paper is the following. In Sec- 
tion [ij we d escribe the first principles method which in 
section |II A is used to calculate the zero field atomic 
structure and Si-H stretch frequencies of the Si(100)- 
H(2xl) surface. The electric fiel d dep endence of the fre- 
quencies is calculated in Section [IB. In Section [II we 



introduce a simple dipole-dipole interaction model for the 
Si-H stretch phonon band, and uses it to find localized 
vibrational states in the presence of electric fields from 
different STM tip geometries. In Section [V we calcu- 



late the lifetimes of the localized states due to incoherent 
lateral diffusion. In section [v| the lifetimes are used to 
model STM induced desorption. 
the results. 



Section VI 



summaries 
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II. STRUCTURE AND VIBRATIONAL 
PROPERTIES OF THE SI(100)-H(2xl) SURFACE. 



In this section we calculate the vibrational frequen- 
cies and the dipole-dipole coupling matrix elements of 
the H vibrations on the Si(100)-H(2x 1) surface. In sub- 
section A we present calculations for the unperturbed 
Si(100)-H(2x 1) surface, and the vibrational and struc- 
tural shifts due to an external planar held are calculated 
in subsection B. 

The hrst principles calculations are based on density 
functional theory [ 19|,^0| within the Generalized Gradient 
Approximation (GGA) |HJ for the exchange-correlation 
energy. Since we only consider filled shell systems, the 
calculations are all non-spin-polarized. Ultra-soft pseudo 
potentials p^] constructed from a scalar-relativistic all- 
electron calculation are used to describe H and silicon(Si) 
p3[ . The wave functions are expanded in a plane- wave 
basis set with a kinetic-energy cutoff of 20 Ry, and with 
this choice absolute energies are converged better than 
0.5 mRy/atom. 

With this approach we find a Si lattice constant of 
5.47(5.43 @)A and bulk modulus of 0.89(0.97 §|) 
Mbar (parentheses show experimental values). For the 
H2 molecule we obtain a bond length of 0.754(0.741 
p4|)A, a bin ding energy (including zero point motion) 
of 4.22(4.52 |24|]) eV and a vibrational frequency of 
4404(4399 |2^])cm _1 . Generally the comparison with ex- 
periment is excellent, and similar theoretical values have 
been found in other studies using the GGA plLp5|. 



A. Zero field properties 



To model the monohydride Si(100) surface at zero field 
we use a (2x1) slab with 12 Si atoms and 6 H atoms. 
The atoms at the bottom surface are bulk like, and their 
dangling bonds are saturated with H atoms. The two 
surfaces are separated by a 7.5 A vacuum region, and 
we use the dipole correction |26f| in order to describe the 
different workfunction of the two surfaces. The surface is 
insulating and we use 2 /c-points in the irreducible part 
of the Brillouin-zone (BZ) for the BZ integrations. Test 
calculations with more dense meshes show that BZ inte- 
gration errors are negligible small. 

Figure 1 shows the atomic structure after relaxation 
of the H atoms and the 4 upper Si layers. Positions of 
the H atoms and the two upper Si layers compare well 
with other studies fl27j-|29| . For the third and fourth layer 
Si atoms we find a small asymmetric relaxation, and to 
our knowledge such relaxations have not been included 
in previous studies. 




FIG. 1. Calculated atomic structure of the Si(100)-H(2x 1) 
surface. 

To obtain the dynamical matrix of the Si-H 
stretch frequency we make H displacements of 
±0.07,0.14, ...,0.35 A in the Si-H bond direction and 
fit a sixth order polynomial to the data points. Since 
the Si-H stretch frequencies are four times higher than 
Si bulk frequencies, the Si substrate acts like a solid 
wall, and we therefore use the H mass Mh when cal- 
culating frequencies from the dynamical matrix. With 
this approach we have calculated the frequency of the 
symmetric stretch lo s and the asymmetric stretch u a at 
three high symmetry points in the surface Brillouin zone 
p?j| . In Table 1 the results are listed together with T 
point frequencies (parenthesis) obtained by infrared spec- 
troscopy (3^] and the comparison between theory and 
experiment is excellent, especially we note that theory 
correctly predicts the splitting, ui s — w a = 11 cm -1 . 

TABLE I. Surface phonon frequencies in cm" 1 for the H 
symmetric stretch ui a and asymmetric stretch <^ a at the F, 
J and J' point in the surface Brillouin zone. Experimental 
r point frequencies are shown in parenthesis. Bracketed J 
and J' point frequencies are obtained from a dipole-dipole in- 
teraction model, with parameters fitted to the two V point 
frequencies (see Section III). 
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We next investigate the bonding potential of the H 
atoms. In Fig. 2 the solid circles show the hydrogen en- 
ergy, £h, for the Si-H bond lengths, z, used in the calcu- 
lation of the dynamical matrix. The data are accurately 
described by a Morse potential 



Eh(z) =E d (e 



— 2a(z — zq ) 



2e 



-a(z-zo) 



), 



(1) 



and from a least-squares fit we obtain the frequency 
luq = 0.26 eV(a = 1.57 A^ 1 ), equilibrium bond length 
z a = 1.50 A, and desorption barrier — 3.4 eV. The 
extrapolated desorption energy coincides with the sur- 
face energy without a H atom plus the energy of a spin- 
polarized H atom. The triangles in Fig. 2 show the total 
energy for large values of z — z$. When the interaction 
between the H atom and the surface becomes weak the 
electrons start to spin-polarize and the data points in 
Fig. 2 show that spin polarization effects become impor- 
tant for z > 2.5 A. 

The inset shows the change in the surface dipole, Ap 
as a function of z (the positive direction is from Si to H). 
The surface dipole increases almost linearly with z and 
the dynamic dipole moment is 7 = dp/dz = 0.6 Debye/A 
(=0.13 e). Modelling the surface dipole by an effective 
charge e* on the H atom and its image charge -e* [^2j , we 
find e* « — 0.07e. The sign of this charge transfer from 
Si to H is in accordance with a higher electronegativity 
of H relative to Si il. 



> 




FIG. 2. Filled circles and open triangles show the change 
in the total energy AEu as a function of the Si-H bond length 
z, and the solid line shows a Morse potential fit to the filled 
circles data points. The horizontal dashed line shows the des- 
orption energy, Ed, as obtained from a spin-polarized H atom 
and the surface without a H atom. To obtain the dynamical 
matrix we fit a sixth order polynomial to the filled circles data 
points. Inset shows the change in the dipole-moment Ap as 
function of z. 



B. Electric field dependent properties 

To model the surface in a planar external field we use 
a (2x 1) slab with 24 Si atoms, 2 H atoms, and a vacuum 
region of 10 A and the external field is modelled using the 
method of Ref . . The Si atoms at the back side of the 
slab are not passivated by H atoms, and dangling bonds 
on these atoms can donate free electrons and holes. In 
this way we take into account the effect of mobile carriers 
p4j| . Other computational details are identical to those 
for the zero field calculation. 




1 2 
E (V/A) 



FIG. 3. The field dependence of (a) the equilibrium Si-H 
bond length Zq, (b) the T point symmetric Si-H stretch fre- 
quency uj a , (c) the r point symmetric-asymmetric splitting 
Au — ujs — ujz- The solid lines show second order polynomials 
obtained by least-squares fit to the data. 

Curves in Fig. 3 show the field dependence of the 
equilibrium Si-H bond length z, the T point symmetric 
Si-H stretch frequency u> s and the T point symmetric- 
asymmetric splitting Alu — lu s — uj a . We first notice that 
all three quantities have an extremum at 1.5 V/A. 
This behaviour can be described by a simple Si-H tight- 
binding model with a field dependent H on-site element 
Jl7j . The extrema occurs at the field where the H and Si 
on-site levels are in resonance, since at resonance the Si-H 
bond is strongest, and therefore the bond length minimal 
and the vibrational frequency maximal. Furthermore, at 
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resonance the H dynamic dipole vanishes, and therefore 
also the part of Aco caused by H-H dipole interactions. 



The three solid lines in Fig. 3 show second order poly- 
nomials obtained by least squares fit to the data. The 
interpolated zero field values of zq and lo s agrees exactly 
with those obtained in section II A , while the interpolated 



Auj zero field value is slightly off. Taking into account 
the quite different slabs used for the two calculations we 
find the agreement fully satisfactory, and note that the 
difference can be taken as a measure of the accuracy of 
the approach. 

Recently, the electric field dependent properties of the 
H/Si(lll)(lxl) surface were calculated by Akpati et al. 
p3| , and they found Stark shifts ~ 30 percent larger 
than in the present calculation, and the extremum in 
bond length and frequency appears for a field of 1 V/A. 
The agreement with the present calculation seems rea- 
sonable, bearing in mind that the Stark shifts are for dif- 
ferent crystallographic surface directions. However, part 
of the difference might be due to the use of a cluster 
geometry and the local spin density approximation in 
Rcf . |33| . The present study is based on a slab geometry 
and the GGA. We expect that the thick slab geometry 
better describes the electric field induced polarization of 
the surface. 
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FIG. 4. The (2x1) surface lattice of H atoms, in units of 
the surface lattice constant aioo = 3.87 A. 



III. STM INDUCED STARK LOCALIZATION 



In this section we will model the collective modes of the 
Si-H stretch vibrations by a set of local oscillators inter- 
acting through dipole forces [ ^5|j36| , and use this model 
to calculate the Stark localization in the external electric 
field from a STM tip. In Fig. 4 is shown the lattice sites 
of the oscillators, corresponding to the positions of the H 
atoms in the (2x1) cell. Each oscillator is described by 
a local frequency uii and a dynamic dipole moment 73. 
The Hamiltonian of the system is given by 



H VJ = hui&u + , _J, 3 (1 " fy), (2) 
r< Ml 



w here r\ is the position of oscillator i, and Xi = 
v / h/2M H uj il i. 



We first consider the zero field case of identical oscil- 
lators with parameters 70, u)q. We find the dispersion 
from a numerical Fourier transform, and at the T point 
we obtain the Hamiltonian 




where Vq — (7q)/(MhWoOioo) an d aioo = 3.87 A is 
the surface lattice constant. The two eigenmodes are 
uj s (T) = u + 9.17Vb, and u} a {T) = uu - 1.{W - Us- 
ing the calculated values of 70 and loq from section 2A, 
we obtain Vb = 0.07 meV and thereby Aut — 5 cm -1 . 
Dipole-dipole interactions can therefore only account for 
half of the dispersion obtained in the the frozen phonon 
calculation. We suggest that the remainder of the split- 
ting is due to a short-range electronic interaction. This 
electronic interaction gives rise to an additional splitting, 
and explains why Auj > at E = 1.5 V/A(see Fig. 3c) 
even though the dipole-dipole interaction vanishes at this 
field. 

To simplify the calculations we will in the following use 
the dipole-dipole interaction modcl(Eq. (||)) to describe 
all the interactions, and determine field dependent pa- 
rameters <x>e and 7e by relating the T-point eigenmodes 
of the model to the calculated frozen phonon values. In 
this way we approximate the effect of the short-range 
electronic interactions by long-range dipole forces. To 
test the accuracy of this approximation we have used the 
model to calculate co s and uj a at the J and J' point in the 
surface Brillouin zone. The result is shown in Table 1 
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and the comparison with the first principles calculation 
is reasonable. 




Surface 



FIG. 5. The geometry used to model the electric field from 
the STM tip, Rt is the tip curvature, s the size of the protru- 
sion on the tip, h the tip-surface distance, and r the radial 
distance from the tip apex. 



We next model the electric field below the STM tip. 
Usually it is found that the tip has a curvature in the 
range 100 A - 1000 A |57],[58| and it is generally accepted 
that the atomic resolution arises from a small protrusion 
or a single atom sticking out of the tip. We use the geom- 
etry in Fig. ^ to model such a tip. Two parameters, the 
tip curvature, Rt, and the protrusion size s, determine 
the tip geometry and we present results for parameters 
in the range, R t = 100-500 A and s = 0-9 A. For the tip- 
sample distance we use, h = 3-7 A, which is the typical 
distance range in STM lithography experiments. To find 
the electric field below the tip the Poisson's equation is 
solved numerically using ANSYS finite element analysis 
p9[ . Curves in Fig. ^ show the radial electric field at 
the surface for a potential difference of 5 V between the 
tip and the surface. Curves in Fig. ^a show the result 
when there is no protrusion on the tip (s = A), and for 
this geometry the electric field attains its half value at 
r rj V Rh. In Fig. ||b results are for a tip with a protru- 
sion of size s, and the protrusion gives rise to a reduced 
electric field below the tip and it decays rapidly around 
the tip apex. The small protrusion changes the electric 
field of the tip very little, and localization of the electric 
field is most pronounced for the large protrusion. Curves 
in Fig. ^|c show the electric field from the geometry with 
Rt — 500 A and s = 6 A at three tip-surface separations 
h = 3,6, and 9 A. The curves show that the field becomes 
more localized when the tip approaches the surface. 
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FIG. 6. Electric field at the surface for the geometry in 
Fig. ^ with a potential difference of 5 V between the tip and 
the surface. We assume that both the tip and the surface are 
metallic. The tip parameters are: (a) s=0 A(no tip protru- 
sion), h=5 A, and R t = 100, 200, 500 A. (b) R t = 500 A, 
h=5 A, and s=0, 3, 6, 9 A. (c) R t = 500 A, s=6 A, and h=3, 
5, 7 A. 

To determine the vibrational states below the tip in 
the presence of the electric field, we set up the Hamil- 
tonian in Eq. (2) for a finite cluster including sites up 
to a cutoff radius r cut and diagonalizes it numerically to 
find the eigenmodes ipa and frequencies u a . There may 
be several localized modes, but we are only interested in 
the localized state with the largest projection p at the 
site directly below the STM tip (r = 0). This state is 
determined using 

p = max[|(^|0)| 2 ], (4) 

a 

where the maximum is over the eigenmodes with 
frequency outside the phonon band, uj a — ujq ^ 
[— 2.8Voj 9.2Vo]- For the spatial electric fields considered 
in this paper the value of p is converged for cluster sizes 
r cut = 50-100 A. 

Curves in Fig. [?] show p and the corresponding vibra- 
tional frequency ui p when the surface is subject to the 
fields of Fig. ||a. The "local E" curve corresponds to the 
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geometry of Persson and Avouris |17| where the electric 
field is localized at r — 0. In this case a localized state 
is split of the phonon band at all negative fields, while 
at positive bias a threshold field of 0.12 V/A is needed 
to obtain localization. For typical fields in H desorption 
experiments, 0.5-1 V/A, the state is completely localized 
at the site below the tip (p — 1) and lo p is similar to the 
frequency, we, of the local oscillator at r — 0. In the 
case of a tip with radius, i? t , a localized state exists for 
nearly all fields, i.e. at positive bias the threshold field 
is 0.03 V/A. The lower positive threshold field compared 
to the "local E" case is obtained because the mode is a 
superposition of several sites with u>i ~ u>e- For typical 
fields in desorption experiments the mode has a substan- 
tial weight, p ~ 0.3, at the site below the tip. 

In Fig. H we show the effect of a small protrusion on 
the STM tip. In this case the spatial localization is im- 
proved, and for fields 0.5-1 V/A we have p ~ 0.8. Thus 
we confirm the results of Persson and Avouris [jl7| , that 
there exists a localized mode in the region below the tip, 
however, it is not completely localized at a single site. 
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FIG. 7. a) The largest projection p at site r = of the 
localized vibrational states when the surface is subject to the 
electric fields of Fig. ^a)(a spherical tip with radius Rt), and 
when the field is localized at r = (solid line). E is the elec- 
tric field at r = 0. b) The vibrational frequency, u p , of the 
localized state relative to the local frequency we at site r — 0. 
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FIG. 8. Similar to Fig. ^a but for the electric fields of (a) 
Fig. |b, and (b) Fig. |c. 



IV. DECAY OF THE LOCALIZED VIBRATION 

Consider an STM experiment where a tunneling elec- 
tron scatters inelastically with the H atom below the tip 
and the H atom is excited into the n = 1 vibrational state 
of the stretch mode. We now consider the decay of such 
an excitation. There are three important time scales, 
the coherent transfer time, r c , the phase relaxation time 
T p h and the energy relaxation time r en . The coherent 
transfer time is the time it takes for the local excitation 
to be transfered into the localized eigenmode below the 
tip, t c ~ h/Aoj w 0.5 ps. Next the eigenmode looses its 
phase due to coupling with a 200 cm" 1 Si phonon pQ 16 
and the phase relaxation time has been measured to be 
Tph ~ 8 ps |^] at room temperature, and t p h s» 75 ps 
|16| at 100 K. Finally the energy of the mode will de- 
cay into the Si substrate via a coupling with three Si-H 
bending modes(600 cm" 1 ) and one 300 cm" 1 Si phonon. 
The time scale for this process is r e „ 10 ns at room 
temperature [ ^p|Jl?jl . 

In the previous section we found a localized eigenmode 
with p 0.8. The excitation at r — will be a super- 
position of this mode and more extended states. After 
t c the extended states have diffused away, thus 80% of 
the initial excitation is in the localized eigenmode, and 
the total probability of finding the initial excitation at 
r = is p 2 ~ 0.6. For t > T p h the excitation can diffuse 
away to the neighbouring H atoms due to dipole-dipole 
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couplings. This is the so-called Forsters mechanism for 
incoherent diffusion, and in the following we will calculate 
the incoherent diffusion rate, w, using Forsters formula 
©0 



w 



|fr w | 9 ^8(wK ? Hdw. 



(5) 



In this equation A® is the spectral function at site i for 
noninteracting H modes (H^ cx dij), but including the 
coupling with substrate phonons which gives rise to the 
phase relaxation. The spectral functions are obtained 
from the noninteracting retarded Greens functions 



-<e(t)<[ci(t),cJ(o)]), 

2ImG°( W ), 



(6) 
(7) 



where cj and Cj are local creation and annihilation oper- 
ators of the stretch mode. The phase relaxation can be 
described approximately by the Hamiltonian J42[ 



(8) 



where rii = b \bi is the projected occupation operator of 
the £1 — 200 cm' 1 Si phonon, and Scj is the change in 
the local frequency when the Si phonon is excited from 
level uton+1. The correlation functions of 7ij have been 
calculated by Persson et al. [|42| 



(n(t)> = n B (fi), 
(n(t)n(0)) = (n)(l + 



(9) 
(10) 



The friction parameter r\ describes the damping of the Si 
phonon, and tib(w) = (e* 3 " — is the Bose occupation 
number and (3 the inverse temperature. 



a) 



b) 
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FIG. 9. Diagrams used for calculating the self energies (a) 



Yi\ and (b) . Solid lines symbolize propagators of the 
stretch mode, while dashed lines are propagators of the Si 
phonon. Wavy lines symbolize the interaction Su>. 



The term gives rise to a small frequency shift, while 
the E^ 2 ) leads to a damping of the mode. Considering 
only the latter term, we find 



r t -(w) 



Ti(u)/h 



( w -^) 2 + r 4 (a>) 2 /4' 

25lu 2 n B (n)[n B (ri) + 1] 
77 (cj - uji) 2 /rj 2 + 1 



(13) 
(14) 



where Ti(cji) is the phase relaxation rate and we have 
used Ti(u)i) 3> rj. Thus the spectral function resembles a 
Lorentzian with Full Width at Half Maximum(FWHM) 
Ti(uji) for cj ~ uiq and it decays as (cj — cjq) 4 in the 
tails. From the experimental dephasing lifetimes |^0jl6 
we obtain Ti(u>i) = l/r p /i ~ 0.7 cm -1 at room temper- 
ature. We estimate the coupling strength using Suji « 
— VLcJi/AEa = —4 cm -1 pij , and the friction parame- 
ter can then be determined from 77 = 2<5cj 2 71b(^)(1 + 
71b(^)) /2Ti(oJi) w 50 cm -1 . The values of 5cj and 77 
obtained in this way are similar to the measured room 
temperature values for Si(lll) p5[. 
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We now use the Matsubara formalism |43|] to obtain 
G° from an perturbation expansion in Sui. We only con- 
sider the two lowest order diagrams shown in Fig. |^, and 
the corresponding self energies are 



e| 1} = H5uj(n), 



,(2) 



= h5uj 2 (n)(l + (n)) 



1 

cj — LOi + irj 



(11) 
(12) 



FIG. 10. Incoherent lateral diffusion rates, w, for the local 
n = 1 vibrational excitation at r = 0, when the surface is 
subject to the electric field of a spherical tip with radius R t 
at h = 5A above the surface(broken curves), or the field is 
localized at r = O(solid). E is the electric field at r = 0.(see 
also Fig. ||) 



7 





a) ; 




n=l 


; ^^^^ /v ? =jUUA, /PjA 

r s=3A 

: S=6A 

r----.-.^..^ - j v=yA 




i i i i i i i i i i i i i i i i i i i - 


; i 1 1 1 1 1 i 1 I i 1 1 1 

r A \ 


.... b) ; 


f ^^<^R =500A, s=6A " v 
' ' /)=3A 

: — /i=6A 

1 /j=9A 


n=l : 




«=2 - 



-1.0 -0.5 0.0 0.5 1.0 

E (V/A) 



of the H atom is well described by a Morse potential. 
The eigenstates of a Morse potential is given by 



ah 1 
(n+ -) 



(16) 



where n takes positive integral values from zero to 
the greatest value for which n + i < \fZM\iEd/ ah. 
For the H potential n = 0,1,... 24 and u{n) = 
0.129,0.378,0.618,0.847,... eV. The anharmonicity is 
substantial and w(2) — w(l) = — w(0) + U, where 
U = —0.010 eV. The frequency of the n = 2 state is out- 
side the phonon band, and this gives rise to a localization 
of the state . The diffusion rate of this state can 

be estimated from Eq. ([l^) by using lu + U for the fre- 
quency at site 0, and the result of such a calculation is 
shown by the three lower curves in figure [□]. The value 
of w is of the same order of magnitude as the room tem- 
perature energy relaxation rate (~ 10 s s _1 ). For n > 2 
the relaxation rate is <C 10 s s _1 . Thus, it is mainly the 
lifetime of the n = 1 excitation which is affected by inco- 
herent diffusion. In the next section we will investigate 
the effect of the reduced lifetime of the n = 1 excitation 
on STM induced desorption. 



FIG. 11. Incoherent lateral diffusion rates, w, for the n = 1 
and n — 2 mode at r = as function of the electric field of 
(a) geometries in Fig. ^p, and (b) geometries in Fig. ^3. 

To obtain the diffusion rate we perform the integration 
in Eq. (||) thus obtaining 



w 



r (wo) + r (wi) 



h 2 f^ rf [^-a;o] 2 + [r K)+r ( Wi )] 2 /4' 



(15) 



In the case where To(u)i) = Tq(luq) this result is similar 
to that of Ref. Q. 

Curves in Fig. show the values of w as obtained 
from Eq. ( fL5| ) when the surface is subject to the same 
electric fields as in Fig. 0. The solid line corresponds to 
the electric field model of Persson and Avouris and sim- 
ilar to Ref. we find w ~ 5 x 10 9 s _1 for typical 
STM fields. The other curves in Fig. |l(] and the curves 
in Fig. [ll] show that for more realistic models of the tip 
electric field the value of w is more than one order of mag- 
nitude larger, and a typical value in an STM experiment 
is w ~ 10 11 s _1 . Thus, the n = 1 vibrational excitation 
at r = will diffuse away very fast to the nearest neigh- 
bour sites in contrast to the result of Persson and Avouris 
||i"8"f . The reason for this is that for a realistic STM geom- 
etry the electric field at r = is not very different from 
the nearest neighbour sites and there is a large diffusion 
rate into these sites. 

For the decay of the n > 1 excitation we have to take 
into account t he an harmonicity of the Si-H bond poten- 
tial. In section II A it was shown that the bond potential 



V. CALCULATION OF THE DESORPTION RATE 

In this section we will calculate the desorption rate, 
R, of the H atom below the STM tip, due to electron 
inelastic scattering through dipole coupling or by reso- 
nance coupling with the Si-H 5c and 6a* resonances. The 
fraction of electrons which scatters inelastically through 
dipole coupling is given by /* p ps (xo/ea ) 2 » 0.001 @. 
The theoretical model we use for calculating the inelastic 
current, /„,_>„, + m , due to resonance coupling has been de- 
scribed in Ref. |ll]|^] . In those works we only considered 
resonance coupling and decay through energy relaxation 



with w c 



l/r eJ 



10" 



and dotted lines in Fig. 



correspond to those results. The dashed lines show the 
result of including dipole coupling and the little differ- 
ence between the dotted and dashed lines justify the ne- 
glect of dipole coupling in our previous studies. The solid 
lines in Fig. 12] show the result of including both dipole 
coupling and lateral diffusion of the n = 1 excitation 
with w = 10 11 s _1 . Defining fx(w) = R(w)/R(w = 0) 
as the suppression of the desorption due to lateral dif- 
fusion of the excitation, we find fi ~ 0.1-0.3 at neg- 
ative bias and /i ~ 0.02-0.08 at positive bias. Using 
w = 10 10 s _1 or w = 10 12 s _1 changes [i less than 10 
percent. This is quite different from the model of Persson 
and Avouris fls| where /i ~ w/w cn . The reason is that 
in our model we include multiple phonon excitations, i.e. 
we use N = 1, 2, 3 in the calculation of the inelastic cur- 
rent Jll]|l^] ■ When the lateral diffusion rate of the n = 1 
level is large, the desorption proceeds via a direct excita- 
tion from n = to n = 2. At negative biases < — 5 V the 
rate of double excitations relative to single excitations is 
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I n ^ n+ 2 / In-tn+i = 0.07-0.15 x (n + 2), while at positive 
biases > 2 V it is 0.015-0.04 x (n + 2). Thus, the larger 
/i at negative bias relative to positive bias is due to a 
higher probability of a multiple excitation. 




Current (nA) 

FIG. 12. Experimental desorption rate Rt as a function of 
tunnel current / for (a) sample bias V& = — 7 V(triangles) 
and —5 V(crosses), (b) sample bias Vb = 2 V(triangles), 
2.5 Vfcircles) and 3 V(squares). Open datapoints are from 
Ref. [ pi|Jl2| and solid data points from Ref. jij. Dashed lines 
show the theoretical results of Ref. [ [ll||l2] | in which inelastic 
dipole scattering and lateral escape of the n = 1 vibrational 
state are not included. Dotted lines show the result when in- 
elastic dipole scattering is included and solid lines show the 
full model where both dipole scattering and a lateral escape 
rate of w — 10 11 s _1 for the n = 1 vibrational state are in- 
cluded. 



VI. SUMMARY 

We have studied the effect of electric field on incoherent 
lateral diffusion of vibrational excitations and its impli- 
cations for STM induced desorption of H from Si(100)- 
H(2xl). We calculated the electric field at the surface 
for realistic STM tip geometries and determined the field 
dependent vibrational properties of the H overlayer based 
on first principles calculations of vibrational Stark shifts 
and dipolc-dipole interaction matrix elements. We found 
that the electric field will localize the vibrational states 
below the STM tip, however, the lifetime of the n = 1 



excitations is short (~ 10 ps) due to incoherent diffusion. 
The diffusion of higher level excitations n > 1 is sup- 
pressed due to anharmonic frequency shifts. The damp- 
ing of the STM induced desorption of H due to the lateral 
escape of the n = 1 excitation depends on the fraction of 
multiple phonon excitation events relative to one phonon 
events in the inelastic scattering process. At low posi- 
tive biases we find a damping of the desorption rate by 
/i ~ 0.02-0.08, while at negative bias \i ~ 0.1-0.3, reflect- 
ing the higher probability of inelastic scattering events 
with a multiple phonon excitation at the negative biases. 
There are no adjustable parameters in our model and the 
calculated desorption rates are in quantitative agreement 
with measured desorption rates. 
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